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The NMR relaxation rates of (quasi-)one-dimensional spin-gapped systems such as the S — 1/2 
spin ladder with a diagonal interaction and the Haldane-gap system are investigated at low temper- 

■ atures in magnetic fields. We show that in the regime where the gap is collapsed by magnetic fields, 
the NMR relaxation rate diverges with decreasing temperature and its exponent shows a behavior 

, characteristic of the model as a function of the magnetic field. We compare our results with recent 

experimental findings for the spin ladder system Cu2(C5Hi2N2)2Ci4 and for the Haldane-gap system 
(CH 3 )4NNi(N0 2 )3. 

hj ; 

■ Low-dimensional quantum spin systems have attracted much attention because of their novel characteristics. Inten- 
sive study has been done both theoretically and experimentally for (quasi-)one-dimensional (ID) spin-gapped systems 
such as the Haldane-gap system, the spin-Peierls system and the spin ladder system. It was shown for the ID 5 = 1/2 

i— h ' spin-gapped systems that the feature of each system appears not in thermodynamic quantities but in the behavior of 
^ , the spin correlation functions when the energy gap is collapsed by external magnetic fields [jij || . Using bosonization 
technique, critical and dynamical properties of the spin correlation functions for the 5=1/2 spin-gapped systems 
in magnetic fields were investigated [ji} |^| and the divergence property of the NMR relaxation rate with respect to 
temperature 1/T\ ~ was discussed in the vicinity of the critical field H Cl , where the gap is collapsed M- By 
combining the numerical diagonalization method with the finite-size-scaling analysis based on conformal field theory, 
the critical properties of several (quasi-) ID spin-gapped systems in magnetic fields were investigated between H C1 and 
the saturation field H C2 , and it was shown that the systems are described by the Tomonaga-Luttinger (TL) liquid 
!> i> i 0> ! !■ Experimentally, the diver gence behavior of 1/Ti of Cu2(CsHi2N2)2Cl4, which can be described as 
the 5=1/2 spin ladder system with a diagonal interaction jlo| O, |l2|], was observed in H Cl < H < H C2 [13]. On 
O \ the basis of the experimental findings, the magnetic-field dependence of the divergent exponent a of \/T\ has been 
discussed quantitatively and compared with some theoretical results (lj |[ [l3| fl4) . 
| In this Letter, we investigate the NMR relaxation rates of (quasi-)lD spin-gapped systems such as the 5=1/2 
. spin ladder system with a diagonal interaction and the Haldane-gap system in magnetic fields H Cl < H < H C2 . Let 
^sD ' us first consider the 5=1/2 spin ladder system with a diagonal interaction in magnetic fields as described by the 
C*~) , following Hamiltonian, 
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where fi(= 1,2) indicates two legs of the ladder, J\ denotes intrachain coupling, J 2 denotes the diagonal interaction 
and H is the applied magnetic field. We set g^B = 1- The interchain coupling Ji is assumed to be positive. Under 
the condition J± ^> | J\\ and | J2I, the model can be mapped onto the ID 5 = 1/2 XX Z model in effective magnetic 
fields |,|l|, 
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n = JY, + SfSf +1 + A5f 5f +1 ] § i> (2) 
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where 5| = 2(5f a - 1) = 2(5? 2 - 1),5+ = -\Z25 S + = V25+ 2 and 5r = -\/25 = V2S^ 2 . The effective 
coupling constants and the effective magnetic field are given by using the original parameters such as J = Ji — j J2, 
AJ = \ Ji + jj 2 and H = H — J± — \,J\ - \ j 2 - When the anisotropy is XY-like, such as — 1 < A < 1 for J > 0, 
the system is gapless and is described by the TL liquid. In this case, the critical fields H Cl and H C2 of the original 
Hamiltonian correspond to the negative and positive saturation fields T-^c, respectively, where H c = | J|(l + A). The 
5 = 1/2 spin ladder with strong interchain coupling (J± > 0) in magnetic fields can be mapped onto the 5 = 1/2 XX Z 
chain with A = 1/2. The 5 = 1/2 bond- alternating spin chain with strong antiferromagnetic coupling can also be 
mapped onto the 5 = 1/2 XXZ chain. Since Cu2(C5Hi 2 N 2 )2Cl4 is believed to be well described by the Hamiltonian 
(|l|) with parameters J\ — 1.0, J 2 = —0.55 and J± = 5.5 fll]] , the magnetic property in H Cl < H < H C2 can be 
investigated by the 5 = 1/2 XXZ chain with the anisotropy A ~ 0.284. 

When the NMR measurement is carried out on the nuclei located at the sites different from the electronic spins, 
relaxation occurs through a dipolar interaction between the nuclear and electronic spins. In this case, the NMR 
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relaxation rate of the ID 5 = 1/2 XX Z model takes the form 

1 If 

— ~ T lim - / dk^sm [x\\ (k,uj) + x± (k, cj)1 , (3) 

where x\\(k,uj) and X_i_(fc, cj) are, respectively, the dynamical longitudinal and transverse spin susceptibilities, and are 
given by 

ikx Aujt 



X||(A,w) = / ds / di(5 2 (x,r)5*(0)) |_ t+0+ e- lte e^, (4) 



Xx(*,w) = / / di(5 + (x,r)^-(0)) | r ^ f+0+ e- ite e*-*. (5) 

In — iJ c < H < H c , the magnetization can be renormalized into one of the phase fields by bosonization treatment. 
Incorporating this phase field in the expressions for spin operators, we obtain the dynamical spin correlation functions 
as 1,0 

(S z (x, r)5 2 (0)) ~ m 2 + C ir - 2 + C 2 r^' cos[tt(1 - 2fh)x], (6) 



(S + (x,t)§-(0)) ~ b 1 r- r >' cos[kx\ + D 2r -^ +1 /^ cos[27rmx], (7) 

with r — \J x 1 + (vt) 2 and v being the velocity of the elementary excitation at T = 0. Note that the asymptotic decay 
of the dynamical spin correlation functions of the original Hamiltonian and the effective Hamiltonian shows the same 
behavior fl7|] . By conformally mapping onto the Matsubara strip as vr±ix = (v/irT)sm[(TrT/v)(vT±ix)], we obtain 
the dynamical spin susceptibility at finite temperature . Paying attention to the temperature dependence of 1 /T\ , 
we obtain \jT x ~ C X T + CfT^'- 1 + Z^T^ -1 + D 2 T'n x + 1 h !C - 1 . Since the universal relation rf = 1 is satisfied in 
the TL liquid [|l9| ^o), the divergence behavior appears only in the second and third terms. Thus, 

^ r ^C 2 T^- 1 +D 1 T r ' x -\ (8) 

where the first term originates in the incommensurate mode of the dynamical longitudinal spin susceptibility and the 
last one originates in the staggered mode of the dynamical transverse spin susceptibility. 

The effective Hamiltonian enables us to calculate the critical exponent rf exactly in — H c < H < H c by means of 
the dressed charge Z(X) as r\ z = 2[Z(A)] 2 |2l|. The dressed charge Z(\) is obtained from the integral equation 
Z{\) = 1 + (l/2vr) d\'K(\ - X')Z(X'), where K(\ - A') = - sin(47)/[sinh(A - A' + i2 7 ) sinh(A - A' - i2j)] for 

27 ^ and K(X - A') = -2/[(A - A') 2 + 1] for 27 = with A = cos(27). The cutoff A is determined by the condition 
for the dressed energy e(±A) = 0, where e(A) is obtained from the integral equation in a given magnetic field as 

e(A) = £ (A) + (l/2vr) j\ d\'K(X — A')e(A') with e (A) = 2H - 2 sin 2 (2 7 )/[sinh(A + i 7 ) sinh(A - i 7 )] for 2 7 ^ and 
£o(A) = 2H — 2/ [A 2 + 1/4] for 27 = 0. Note that the magnetization is obtained by the use of the dressed charge as 
m = 1/2-(1/2tt) j\ dXZ(X)q a (X) with q Q (X) = sin(27)/[sinh(A + i 7 ) sinh(A-i 7 )] for 2 7 ^ and q (X) = l/[A 2 + l/4] 
for 2 7 = 0. In the case of < A < 1, the exponent rf- thus obtained varies convexly downward as a function of 
the magnetic field and 1 < rf < 2, while in the case of —1 < A < 0, rf varies convexly upward and rf > 2, with 
rf being 2 at H = ±H C in both cases. The behavior of rf can be obtained from the universal relation rf rf = 1. 
The results for rf and rf have also been confirmed by the numerical results for several systems such as the S = 1/2 
bond-alternating chain ||, the S = 1/2 spin ladder M and the S = 1/2 spin ladder with a diagonal interaction ||. 

On the basis of these observations, we conclude that the divergence property of the NMR relaxation rate 1/T\ of the 
system, which can be mapped onto the S = 1/2 XX Z chain with an XF-like anisotropy, is caused by the staggered 
mode of the dynamical transverse spin susceptibility 

— - DxT^- 1 = D x T- a . (9) 

The results for the exponent a obtained from the effective Hamiltonian and obtained numerically from the original 
Hamiltonian S are shown in Fig. 1 for the 5 = 1/2 spin ladder system and the 5 = 1/2 spin ladder with a diagonal 
interaction corresponding to (^(CsH^^^CLi. The parameters for the numerical calculation of the former system 
are J\ = 1.0, J2 = and J j_ = 5.0, and those for the latter system are indicated previously ||. The magnetization of 
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FIG. 1: Exponents of the NMR relaxation rate 1/Tj ~ D 1 T~ a for the 5=1/2 spin ladder system and the S = 1/2 spin 
ladder with a diagonal interaction corresponding to Cu2(C5Hi2N2)2CU are shown as a function of magnetization m. The solid 
line and the broken line represent the results from the ID S = 1/2 XX Z model with A = 1/2 and A = 0.284, respectively. 
Diamonds and crosses are the results from the numerical data of ref.8. 



the original Hamiltonian m and that of the effective Hamiltonian rh satisfy the relation 2m = m + 1/2. The results 
by both methods show qualitatively the same behavior. In particular, the agreement becomes noticeable in the larger 
m region. In the magnetic field near H C2 , the system is nearly fully polarized and the many-body effect may be 
suppressed, leading to the good agreement between the results of both methods. In the magnetic field near H Cl , on 
the other hand, a quantum fluctuation which relates to a gap may still survive, leading to the difference between the 
results by both methods. 

The result for the S = 1/2 spin ladder system is different from that obtained based on the bosonization technique, 
where the exponent a = 1/2 at m = and does not change with m in the vicinity of m = Q. Furthermore, the 
effective TL parameter at H = H Cl is K = 1/2 in ref. 1, whereas K — r) z /2 = 1 in our result. At H = H Cl and H C2 , 
the system is in the empty or fully occupied band in the fermionic language, which leads K = 1. Thus we conclude 
that in the 5=1/2 spin ladder system, K and a vary convexly downward as a function of the magnetic field with 
K = 1 and a — 1/2 at H — H Cl and H C21 irrespective of the strength of interchain coupling, although the magnetic 
field where K and a take minimum values may shift in the case of weak interchain coupling. 

We discuss further the magnetic-field dependence of 1/Ti in other systems such as the Haldane-gap system and the 
S = 1/2 bond-alternating spin chain with a second-neighbor interaction. The latter system can be described by the 
Hamiltonian ([!]). With the parameters J± = 2.06, J\ = 0.48 and J<i = 1.94, the Hamiltonian (|l|) can be an effective 
model to describe the magnetic properties of CuGeC>3 which exhibits the spin-Peierls transition p2| . 

These two systems seem to be difficult to map onto the effective S — 1/2 XX Z chain. The critical properties of 
these systems in H c ,< H < H C2 were investigated numerically, and it was revealed that the systems are described 
by the TL liquid f|, ||, ■ The ID spin systems described by the TL liquid can be mapped onto the interacting boson 
system, where the density operator of the boson corresponds to S* of the effective ID S = 1/2 XX Z model and the 
creation and annihilation operators of the boson correspond to S . The dynamical correlation functions of these 
boson operators take the same asymptotic forms as (|6|) and (]?]). Therefore, the result for 1/Ti shown in (||), which 
has been derived based on (^J) and (Q), holds generally in ID spin systems described by the TL liquid. Then, we 
investigate the divergence behavior of 1/Ti in ID spin systems described by the TL liquid, which cannot be mapped 
onto the effective ID S = 1/2 XXZ model, using the result (||). 

We first consider the Haldane-gap system. The numerical results indicate that rf < 1/2 and rf > 2 with rj x = 1/2 
and rf = 2 at H Cl and H C2 [Q, ||. Judging from these results and (||), we conclude that the divergence behavior of 
1/Ti is caused by the staggered mode of the dynamical transverse spin susceptibility and 1/Ti is expressed as shown 
(|) (23). The results for a in the Haldane-gap system are shown in Fig. 2 using the numerical results in ref. 4. 



The exponent a increases from 1/2, and takes a maximum of ~ 0.65 around m ~ 0.4. At H = H C2 , a — 1/2. Quite 
recently, the divergence behavior of 1/Ti in the Haldane-gap system (CH3)4NNi(N02)3 was observed experimentally 
at low temperatures ^24j| . The exponent of 1/Ti with respect to temperature has been estimated as 0.5 at m ~ and 
as 0.7 at m ~ 0.2. These experimental results agree with the results shown in Fig. 2 quantitatively. 

We next consider the S — 1/2 bond-alternating spin chain with a second- neighbor interaction with the parameters 
Jj_ = 2.06, J\ = 0.48 and J2 = 1.94. The critical exponents rf and rf show that rf > 1 and rf < 1 in 0.1 < m < m c , 
and \ j2 < rf < 1 and 1 < rf < 2 in m c < m < 1/2, where m c ~ 0.35 and rf — rf — 1 at m — m c 0. From 
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FIG. 2: Exponent of 1 jT\ ~ Z)iT a for the Haldane-gap system is shown as a function of m. The inset shows the H dependence 
of a. 

(§), we conclude that 1/T X - C , 2 T - ( 1 - 7 '*) in 0.1 < m < m c and l/T x ~ D 1 T~( 1_T? *) in m c < m < 1/2. The results 
for the divergent exponent a in this system are shown in Fig. 3 using the numerical results in ref. 7. The results 
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FIG. 3: Exponent of 1/Ti for the S =1/2 bond-alternating spin chain with a second-neighbor interaction is shown as a function 
of m. The parameters are J± = 2.06, Ji = 0.48 and Ji = 1.94. The inset shows the H dependence of the exponent. 

indicate that the leading contribution of the spin correlation function is the incommensurate mode of (S z (x, t)S z (0)) 
in 0.1 < to < m c , and is the staggered mode of (S + (x, t)S~ (0)) in m c < m < 1/2. At H = H Cl , the system is 
in the empty or fully occupied band and K = t] z /2 = 1. Therefore, it is expected that near H = H Cl , the leading 
contribution is the staggered mode of (S + (x, t)S~ (0)) and that the dominant mode changes again at the magnetic 
field between H Cl and 0.1. In the region where the longitudinal incommensurate mode is dominant, the system has a 
tendency to form an incommensurate spin order with period 1/m. This tendency may relate to the property of the 
incommensurate phase of CuGeC>3 in magnetic fields |2J| . 

The behavior of r] x and rf in 0.1 < m < m c is most likely brought about by the long-range interaction: the 
second-neighbor antiferromagnetic interaction. Hence, in systems such as the S =1/2 n-leg spin ladder system with 
diagonal interactions, the incommensurate mode of the longitudinal spin correlation function can be dominant in an 
intermediate range between H Cl and H C2 , while the staggered mode of the transverse one can be dominant near H Cl 
and H C2 . This characteristic behavior is detectable by the NMR relaxation rate in magnetic fields, if the probed 
nucleus is at a site different from electronic spins. 

Finally, we comment on the NMR relaxation rate of the Heisenberg alternating-spin chain with S = 1 and S =1/2 
in magnetic fields. It was shown numerically that the system belongs to the universality class of the TL liquid when 
1/4 < to < 3/4 p6|. Using the numerical results for if and rf |26| and (ph, we see that the divergence behavior of 
1 /T\ is dominated by the staggered mode of the dynamical transverse spin susceptibility, and the exponent seems to 
take a minimum of ~ 0.4 around to = 0.4. 
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In summary, we have systematically investigated the NMR relaxation rates of (quasi-) ID spin-gapped systems at 
low temperatures in H C1 < H < H C2 . The results are summarized in Table I. The divergent exponent of 1 /Ti with 
respect to temperature shows a feature characteristic of the model, reflecting the nature of the TL liquid. 

We would like to thank H. Frahm, T. Goto, N. Kawakami and A. Tomiyama for their useful comments and valuable 
discussion. This work was supported by Grant-in-Aid no. 10640344 for Scientific Research from the Ministry of 
Education, Science, Sports and Culture, Japan. 



TABLE I: The divergent exponents of 1/Ti in magnetic fields are summarized. AF (F) Ladder denotes the 5=1/2 lad- 
der with antiferromagnetic (ferromagnetic) interchain couplings, AF-AF(-F) BA denotes the S = 1/2 antiferromagnetic- 
antiferromagnetic(-ferromagnetic) bond-alternating chain and Frustrated BA denotes the system with J± = 2.06, Ji = 0.48 
and J 2 = 1.94. 



Model 


at H C1 and H C2 


Dominant mode 


H C1 < H < H C2 


Cu 2 (C5Hl2N 2 )2Cl4 


1/2 


transverse stag 


gered mode 


convex downward 


AF Ladder 


1/2 


transverse stag 


gered mode 


convex downward 


AF-AF BA 


1/2 


transverse stag 


gered mode 


convex downward 


Alternating-spin chain 


1/2 


transverse stag 


gered mode 


convex downward 


Haldane Gap 


1/2 


transverse stag 


gered mode 


convex upward 


AF-F BA 


1/2 


transverse stag 


gered mode 


convex upward 


F Ladder 


1/2 


transverse stag 


gered mode 


convex upward 


Frustrated BA 


1/2 


dominant mode changes 


depending on the mode 
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